A cost‐effective curvature calculation approach for interfacial flows on unstructured meshes

We present a simple and cost‐effective curvature calculation approach for simulations of interfacial flows on structured and unstructured grids. The interface is defined using volume fractions, and the interface curvature is obtained as a function of the gradients of volume fractions. The gradient computation is based on a recently proposed gradient recovery method that mimicks the least squares approach without the need to solve a system of equations and is quite easy to implement on arbitrary polygonal meshes. The resulting interface curvature is used in a continuum surface force formulation within the framework of a well‐balanced finite‐volume algorithm to simulate multiphase flows dominated by surface tension. We show that the proposed curvature calculation is at least as accurate as some of the existing approaches on unstructured meshes while being straightforward to implement on any mesh topology. Numerical investigations also show that spurious currents in stationary problems that are dependent on the curvature calculation methodology are also acceptably low using the proposed approach. Studies on capillary waves and rising bubbles in viscous flows lend credence to the ability of the proposed method as an inexpensive, robust, and reasonably accurate approach for curvature calculation and numerical simulation of multiphase flows.     

风力机叶片气动噪声的影响参数*

随着风力机的大型化，风电机组对环境的影响不容忽视，必须对风力机气动噪声进行预测和控制。选取基于NACA、DU翼型的某风力机叶片作为研究基准，采用修正BPM半经验模型计算叶片的气动噪声特性，通过改变翼型族、弦长、机组运行状态、风切变指数、来流风向参数，研究叶片外形几何参数、机组运行工况对叶片气动噪声源的影响。计算结果从多个角度总结出水平轴风力机叶片气动噪声的变化规律，为开发高效低噪风电叶片提供参考。     

Comparison of implicit and explicit AUSM‐family schemes for compressible multiphase flows

This paper makes the first attempt of extending implicit AUSM‐family schemes to multiphase flow simulations. Water faucet, air–water shock tube and oscillating manometer problems are used as benchmark tests with the generic four‐equation two‐fluid model. For solving the equations implicitly, Newton's method along with a sparse matrix solver (UMFPACK solver) is employed, and the numerical Jacobian matrix is calculated. Comparison between implicit and explicit AUSM‐family schemes is presented, indicating that similarly accurate results are obtained with both schemes. Furthermore, the water faucet problem is solved using both staggered and collocated grids. This investigation helps integrate high‐resolution schemes into staggered‐grid‐based computational algorithms. The influence of the interface pressure correction on the simulation results is also examined. Results show that the interfacial pressure correction introduces numerical dissipation. However, this dissipation cannot eliminate the overshoots because of the incompatibility of numerical discretization of the conservative and non‐conservative terms in the governing equations. The comparison of CPU time between implicit and explicit schemes is also studied, indicating that the implicit scheme is capable of improving the computational efficiency over its explicit counterpart. Copyright © 2014 John Wiley & Sons, Ltd.     

适用于等熵流动的交错拉氏Godunov方法

为消除传统单元中心型Godunov方法在求解稀疏波问题时的非物理过热现象，发展一种适用于等熵流动的交错拉氏Godunov方法.主要的特征是采用速度与热力学变量交错分布的形式，避免在单元内进行速度平均，从而消除由于动量平均过程导致的动能耗散.与传统的von Neumann型交错网格方法相比，网格的边界通量由节点处的多维黎曼求解器提供，克服了多维人工粘性选取带来的困难.为减少多维黎曼求解器在求解稀疏波问题时的非物理熵增，给出稀疏波出现的合理判据，从而保证了热力学关系式的满足.数值实验表明：该方法能很好地消除稀疏波的过热现象，同时在求解激波问题时又能保持与传统单元中心型拉氏方法相同的激波捕捉能力.     

New Sixth-Order Compact Schemes for Poisson/Helmholtz Equations

Some new sixth-order compact finite difference schemes for Poisson/Helmholtz equations on rectangular domains in both two- and three-dimensions are developed and analyzed. Different from a few sixth-order compact finite difference schemes in the literature, the finite difference and weight coefficients of the new methods have analytic simple expressions. One of the new ideas is to use a weighted combination of the source term at staggered grid points which is important for grid points near the boundary and avoids partial derivatives of the source term. Furthermore, the new compact schemes are exact for 2D and 3D Poisson equations if the solution is a polynomial less than or equal to 6. The coefficient matrices of the new schemes are $M$-matrices for Helmholtz equations with wave number $K≤0,$ which guarantee the discrete maximum principle and lead to the convergence of the new sixth-order compact schemes. Numerical examples in both 2D and 3D are presented to verify the effectiveness of the proposed schemes.     

航空发动机压气机叶片振动频率与温度的关系

以某型航空发动机压气机第六级叶片为例，研究了温度对叶片振动频率的影响。发现由于温度的升高导致了材料弹性模量的下降，而弹性模量的下降又导致了叶片在发动机一定的转速范围内，除一阶振动外，其余各阶振动的动频均小于静频，因此与通常的工程处理方法不同，研究压气机叶片的动频时，应考虑温度的影响。     

风机叶片包容性分析

应用有限元方法对丢失的失效风机叶片撞击机匣的非线性瞬态响应进行了数值计算研究,模拟了叶片的包容过程,分析了包容过程中叶片、机匣的变形与能量变化。结果表明:有限元方法能较好的模拟风机叶片丢失后撞击机匣的过程,该型风机机匣对叶片具有包容性。研究结果对于进行风机叶片包容性设计具有重要参考价值。     

基于SIMPLER算法的多重网格方法研究

以二维方腔顶盖驱动流为模型,将多重网格方法和SIMPLER算法进行耦合,对不同雷诺数下多重网格加速SIMPLER算法和SIMPLER算法的计算效率进行了对比,数值计算表明:多重网格加速SIMPLER算法不仅能够解决SIMPLER算法不能准确模拟较高雷诺数流场的问题,而且其计算效率远远高于SIMPLER算法.本文也对松弛因子的选取、多重网格实现形式以及网格层数对多重网格加速SIMPLER算法的影响进行了研究,从而为多重网格加速SIMPLER算法的实施提供了计算技术.     

涡轮叶片前缘阵列冲击冷却流动及传热特性数值研究

本文采用SST湍流模型模拟了类前缘通道内蒸汽射流阵列冲击冷却的流动与传热特性,分析了雷诺数(Re=10000~50000)、孔径比(d/H=0.5~0.9)和孔间距比(S/H=2~6)对流动及传热性能的影响规律,得到了相应的传热和摩擦关联式。结果表明:在不同雷诺数下,d/H从0.5到0.9变化时,通道压力损失系数降低了76%~79%,靶面平均努塞尔数降低了45%~49%;S/H从2增至6时,通道压力损失系数增加了1.64~1.92倍,靶面平均努塞尔数增加了54%~64%;增大d/H、减小S/H可有效提高类前缘通道蒸汽冲击冷却的综合热力系数。本文研究结果可为未来先进燃气轮机高温涡轮叶片蒸汽冷却结构的设计提供参考和借鉴。     

Lattice Boltzmann simulation of two-sided lid-driven flow in a staggered cavity

In this article, the lattice Boltzmann method is employed in order to explore incompressible fluid flow inside a two-sided lid-driven staggered cavity. Results of the lattice Boltzmann simulation for antiparallel motion of lids are compared with the data from existing literature. For parallel motion of lids, the characteristics of flow pattern for a variety of Re numbers (50–3200) are presented. An asymmetric steady-state flow pattern for parallel motion of lids is obtained.